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We  investigate  experimentally  and  theoretically  the  be¬ 
havior  of  Coulomb  blockade  (CB)  peaks  in  a  magnetic  field 
that  couples  principally  to  the  ground-state  spin  (rather  than 
the  orbital  moment)  of  a  chaotic  quantum  dot.  In  the  first 
part,  we  discuss  numerically  observed  features  in  the  magnetic 
field  dependence  of  CB  peak  and  spacings  that  unambiguously 
identify  changes  in  spin  S  of  each  ground  state  for  successive 
numbers  of  electrons  on  the  dot,  N.  We  next  evaluate  the 
probability  that  the  ground  state  of  the  dot  has  a  particu¬ 
lar  spin  S,  as  a  function  of  the  exchange  strength,  J,  and 
external  magnetic  field,  B.  In  the  second  part,  we  describe 
recent  experiments  on  gate-defined  GaAs  quantum  dots  in 
which  Coulomb  peak  motion  and  spacing  are  measured  as  a 
function  of  in-plane  magnetic  field,  allowing  changes  in  spin 
between  N  and  N  +  1  electron  ground  states  to  be  inferred. 

I.  INTRODUCTION 

In  the  absence  of  interactions,  electrons  populate  the 
orbital  states  of  a  quantum  dot  or  metallic  grain  in  an 
alternating  sequence  of  spin-up  and  spin-down  states,  in 
accordance  with  the  Pauli  picture.  In  this  case,  the  total 
spin  of  the  dot  is  either  zero,  when  the  number  of  elec¬ 
trons  N  is  even,  or  one  half  at  odd  N.  No  higher  spin 
states  appear.  It  is  well  known  that  in  the  presence  of 
interactions,  this  simple  scheme  can  be  violated.  Such 
deviations  from  simple  even-odd  filling  may  appear  as 
Hund’s  rules  in  atomic  physics,  or  be  as  extreme  as  com¬ 
plete  spin  polarization  leading  to  ferromagnetism  by  the 
mechanism  of  the  Stoner  instability. 

By  measuring  electron  transport  through  a  weakly  cou¬ 
pled  quantum  dot  at  low  temperature  and  bias  voltage 
(compared  to  the  quantum  level  spacing)  one  can  effec¬ 
tively  study  differences  in  ground  state  (GS)  properties  of 
the  dot  at  successive  electron  number,  providing  a  means 
of  investigating  the  actual  filling  scheme  jjj.  A  stan¬ 
dard  experimental  approach  [jlj  is  to  measure  conduction 
through  the  dot  via  two  tunneling  leads  as  a  function  of 
a  voltage,  Vg ,  applied  between  the  dot  and  a  third  elec¬ 
trically  isolated  electrode,  known  as  a  gate.  Because  the 
dot  GS  energy,  En,  at  a  fixed  electron  number  is  influ¬ 
enced  by  Vg,  the  gate  can  be  used  to  set  the  number  of 
electrons  on  the  dot.  The  fact  that  a  large  Coulomb  en¬ 
ergy  is  needed  to  add  a  single  electron  to  the  dot  typically 
suppresses  conduction  through  the  dot  in  the  tunneling 


regime,  where  whole  charges  must  tunnel  for  conduction 
to  occur.  This  effect  is  known  as  the  Coulomb  blockade 
(CB).  However,  at  specific  values  of  gate  voltage,  denoted 
VgN\  where  the  condition  En- i{Vg)  =  Ejy{Vg)  is  satis¬ 
fied,  conductance  increases  dramatically,  clearly  marking 
this  degeneracy  condition  in  an  experimental  trace.  The 
position,  VgN\  of  the  Nth  peak  in  the  conductance  is  pro¬ 
portional  to  =  En( 0)  —  -Ejv-i(O).  Accordingly,  the 
distance  between  successive  peaks  is  =  hn  —  hn- i- 
(We  take  the  constant  of  proportionality  converting  gate 
voltage  to  dot  energy  to  be  unity  for  the  theoretical  dis¬ 
cussion;  experimentally,  this  constant  can  be  readily  mea¬ 
sured,  for  instance,  by  comparing  the  influence  of  bias 
and  gate  voltages.) 

In  the  last  few  years  the  GS  spin  of  a  variety  of  nanos¬ 
tructures,  including  metallic  grains  semiconducting 
quantum  dots  |(|,  and  carbon  nanotubes  |7|,[|,  have 
been  investigated  using  CB  peak  motion  in  a  magnetic 
field.  If  one  neglects  the  magnetic  field  coupling  to  the 
orbital  degrees  of  freedom,  one  expects  the  field  to  mani¬ 
fest  itself  only  through  the  Zeeman  splitting,  resulting 
in  a  shift  of  the  GS  energy  by  gfisSB,  where  is 
the  Bohr  magneton  and  S  denotes  the  GS  total  spin. 
For  2D  quantum  dots,  orbital  coupling  can  be  strongly 
suppressed — though  not  eliminated  entirely — by  orient¬ 
ing  the  field  strictly  in  the  2D  plane.  This  is  the  exper¬ 
imental  approach  that  will  be  described  in  Section  III 
of  this  paper.  On  the  other  hand,  for  ultrasmall  grains 
and  nanotubes,  it  is  reasonable  to  ignore  orbital  coupling 
for  any  field  direction  since  for  practical  magnetic  fields 
the  total  flux  cutting  the  structure  is  much  less  than  the 
quantum  of  flux. 

For  semiconducting  quantum  dots  |i|  the  Zeeman 
splitting  for  S'  =  1/2  becomes  comparable  to  the  mean 
single  electron  level  spacing 

=  {si) ;  Si  =  Si  —  £j+i,  (1) 

(here  et  denotes  the  orbital  energy  of  the  one-electron 
orbital  state  i,  and  (. . .)  stands  for  an  average  over  dif¬ 
ferent  levels)  at  B  ~  1 T  for  the  GaAs  dots  is  Refs.  [[i|[| 
and  B  ~  10  T  for  the  Si  dot  in  Ref.  || .  At  lower  fields, 
the  splitting  seldom  exceeds  s, ,  and  one  might  expect  a 
simple  even/odd  Pauli  picture.  In  this  case,  CB  peak  po¬ 
sitions  would  move  with  B  as  consecutive  pairs,  creating 
a  pattern  of  alternating  downward  and  upward  moving 
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peaks,  as  seen  in  Fig.  |l|a.  The  trajectories  of  the  pair  of 
peaks  would  form  straight  lines  with  the  same  slope  of 
ggb/2.  The  peak  spacings,  A^,  will  have  a  similar  pat¬ 
tern  with  slopes  of  gg Spin-orbit  interaction  can  lead 
to  fluctuations  in  the  g  factor,  resulting  in  fluctuations 
in  the  slopes  of  the  lines,  [^],^0|,  but  will  not  change  the 
pattern  of  downward  and  upward  moving  peaks.  Once 
ggbB  exceeds  a  particular  spacing  si:  the  peaks  should 
cross. 

Experimentally,  the  Pauli  picture  appears  to  work  well 
for  metallic  grains  |Q,  while  for  carbon  nanotubes  the 
situation  is  less  clear:  one  published  experiment  finds 
simple  even/odd  filling  [Q,  and  one  finds  a  more  com¬ 
plicated  scheme  [|j.  For  multielectron  semiconducting 
dots  the  data  are  somewhat  confusing,  but  ap¬ 

pear  to  disagree  with  the  simple  even/odd  filling  scheme 
presumably  as  a  result  of  electron-electron  in¬ 
teractions.  In  symmetric,  few-electron  vertical  dots  l[i~2| 
one  may  again  recover  relatively  simple  behavior  for  the 
first  few  electrons,  with  a  well-understood  appearance 
of  Hund’s  rules.  These  conclusions  are  based  on  exper¬ 
iments  where  orbital  magnetism  dominates  spin.  How¬ 
ever,  for  few  electrons  (i.e.,  N  <  10)  the  filling  scheme 
can  be  readily  interpreted,  allowing  the  full  spin  structure 
to  be  mapped  out  as  a  function  of  N  and  B  ]l2[ . 

Figure  [j]  shows  CB  peak  motion  and  spacing  as  a  func¬ 
tion  of  magnetic  field  (coupling  only  to  spin)  for  both  a 
noninteracting  electron  model  and  an  interacting  elec¬ 
tron  model  presented  in  detail  in  Section  II.  The  behav¬ 
ior  of  this  model  appears  more  similar  to  the  experimen¬ 
tally  observed  behavior  in  semiconducting  quantum  dots 
(see  Fig.  |J.  In  some  cases  peak  positions  (g n )  move  as 
function  of  the  magnetic  field  in  bunches  of  two  in  the 
same  direction  and  with  the  same  slope,  while  their  cor¬ 
responding  spacings  (A^)  do  not  change  as  function  of 
the  magnetic  field.  Moreover,  not  every  change  in  the 
direction  of  the  motion  of  a  peak  can  be  explained  as  a 
crossing  of  two  orbital  levels.  Thus,  a  description  of  the 
system  beyond  the  Pauli  picture  is  needed. 

We  will  present  a  theory  that  goes  beyond  the  simple 
Pauli  picture  and  explains  the  puzzling  features  of  the 
magnetic  field  dependence  of  the  conductance  peaks  po¬ 
sitions  as  a  result  of  spontaneous  spin  polarization  of  the 
electrons  in  the  dot  0-  The  possibility  of  spontaneous 
magnetization  in  quantum  dots  has  been  considered  in 
different  contexts  [  i~i|  fiojj .  For  instance,  it  has  also  been 
suggested  as  an  explanation  |n|  for  the  absence  of  a  bi- 
modal  distribution  in  the  conductance  peak  spacings  [|l] . 
Such  a  magnetization  was  also  predicted  to  lead  to  kinks 
in  the  parametric  motion  of  the  peaks,  (e.g.,  due  to  or¬ 
bital  effects  of  a  perpendicular  magnetic  field)  El¬ 
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FIG.  1.  Typical  conductance  dependence  on  the  gate  volt¬ 
age  and  magnetic  field,  B,  with  and  without  exchange  in¬ 
teraction:  (a)  J  =  0,  (b)  J  =  0. l<5i .  g'N  is  determined  as 
gN  =  gN  —  e2 /C  and  taken  in  arbitrary  units.  A^  for  each 
N  is  shifted  arbitrarly  by  Nx  constant  to  avoid  overlapping. 

II.  THEORY 

We  begin  the  theoretical  discussion  by  reviewing 
the  considerations  determining  the  general  form  of  the 
Hamiltonian  describing  the  properties  of  a  chaotic  dot  in 
the  metallic  regime.  A  more  complete  description  may 
be  found  in  Ref.  0-  The  single  electron  spectrum  £j 
is  characterized  by  the  mean  level  spacing  <5i  and  the 
Thouless  energy  Et  ~  Ti/terg,  where  terg  is  the  time  it 
takes  for  a  classical  electron  to  cover  the  energy  shell  in 
the  single-particle  phase  space.  For  diffusive  and  ballis¬ 
tic  systems  terg  equals  to  L2  /  D  and  L/vf  respectively, 
where  vf  is  the  electronic  Fermi  velocity,  L  is  the  system 
size,  and  D  is  the  diffusion  coefficient.  The  dimension¬ 
less  conductance,  gr  =  Et/Si,  in  the  metallic  regime 
is  large,  i.e.,  gr  ^  1-  In  this  regime  the  statistics  of 
the  single  electron  spectrum  on  scales  smaller  than  Et 
are  well  described  by  random  matrix  theory  (RMT)  [p2|, 
which  gives  a  quantitative  description  of  the  phenomenon 
of  level  repulsion.  For  an  ensemble  of  M  x  M  matrices 
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(M  — >  oo )  with  random  and  independent  elements,  the 
probability  density  of  a  realization  of  the  spectrum  £;  is 
given  by  p§: 


P(£i)  oc  exp 


(2) 


where  (3  is  equal  to  1,  2  or  4  for  the  orthogonal,  unitary 
and  symplectic  ensembles  respectively.  The  orthogonal 
(unitary)  RM  ensemble  corresponds  to  weakly  disordered 
dots  with  preserved  (violated)  time  reversal  symmetry 
and  negligible  spin-orbit.  We  assume  that  the  dot  is  small 
enough  for  the  spin-orbit  interaction  to  be  negligible  and 
thus  avoid  discussion  of  the  symplectic  ensemble  (/3  =  4) 


The  second  step  is  to  consider  effects  of  electron- 
electron  interactions.  For  the  simplest  case  of  short-range 
interactions,  we  write 


Hint(r)  =  A  S,V6(r),  (3) 

where  V  ot  Ld  is  the  volume  of  the  dot  (in  d  dimensions) 
and  A  is  a  dimensionless  coupling  constant  characteriz¬ 
ing  the  strength  of  the  interaction.  Matrix  elements  of 
this  interaction,  in  the  basis  of  eigenstates  g>i(r)  of  the 
noninteracting  Hamiltonian,  are  given  by 

Ml\  =  \5xv  J  drg}*{r)ip*{r)ipk{f)gn{r).  (4) 

It  is  important  to  note  that  the  statistical  properties  of 
the  interaction  matrix  elements  M ^  are  completely  de¬ 
termined  by  the  statistical  properties  of  the  single  elec¬ 
tron  eigenstates  tpi{r)  and  cannot  be  chosen  in  an  arbi¬ 
trary  fashion  as  in  the  random  interaction  model  [po|. 
In  fact,  in  the  limit  gx  — >  oo  the  diagonal  matrix  el¬ 
ements  (*,  j,  k,  l  pairwise  equal)  are  self-averaged  by  the 
space  integration,  Eq.  (||),  and  thus  show  no  level-to-level 
or  sample-to-sample  fluctuations  0-  In  the  same  limit 
the  off-diagonal  matrix  elements  Eq.  (^[)  turn  out  to  be 
negligable. 

Using  the  statistical  properties  of  the  interaction  ma¬ 
trix  elements  it  turns  out  that  a  large  class  of  disordered 
metallic  dots  can,  under  very  general  conditions,  be  de¬ 
scribed  by  a  remarkably  simple  Hamiltonian  with  only 
three  coupling  constants,  which  do  not  fluctuate.  In  the 
limit  of  large  gr ,  the  interaction  part  of  the  Hamiltonian 
corresponds  to 

Hint  =  ECN 2  -  J(S)2  +  A bcsT^T,  (5) 

(terms  linear  in  N  are  allowed,  but  they  can  be  included 
into  the  one-particle  part  of  the  Hamiltonian)  where 
N  is  the  number  operator,  S  is  the  spin  operator  and 
T  =  (ci|  annihilates  an  electron  in  the  ith  sin¬ 

gle  electron  orbital  with  spin  f)-  In  the  simple  model 
with  short-range  interaction  and  preserved  time  reversal 


(/3  =  1)  the  above  coupling  constants  have  the  following 
form: 


Ec  =  J  =  2A<5i;  A  Bcs  =  A<5i .  (6) 

If  time  invariance  is  broken  {(3  =  2),  A  bcs  =  0  since  the 
operator  T  is  incompatible  with  the  symmetry.  Note  that 
the  interaction  Hamiltonian  Eq.  (||)  represents  a  particu¬ 
lar  model,  and  the  expressions  for  the  coupling  constants 
Eq.  (jcij)  are  valid  only  for  this  model.  At  the  same  time, 
the  effective  interaction  Hamiltonian,  Eq.  (||),  is  more  ro¬ 
bust  and  depends  only  on  the  symmetries  of  the  problem 
(for  instance,  on  the  absence  of  the  spin-orbit  scattering) 
and  on  the  condition  gr  1- 

The  first  two  terms  in  Eq.  (|j)  represent  the  dependence 
of  the  energy  of  the  dot  on  the  total  number  of  the  elec¬ 
trons  and  on  the  total  spin  respectively.  They  commute 
with  each  other  and  with  the  single-particle  part  of  the 
Hamiltonian.  Therefore,  all  states  of  the  grain  can  be 
classified  by  N  and  S.  The  term  proportional  to  A  bcs 
appears  only  in  the  orthogonal  case  (/3  =  1).  Provided 
that  A  bcs  <  0  this  term  leads  to  the  superconducting 
instability.  Superconducting  correlations  are  suppressed 
by  the  magnetic  field,  and  thus  do  not  exist  for  (3  =  2. 

Thus,  the  general  form  of  the  Hamiltonian  describing 
electrons  in  a  chaotic  dot  which  is  not  superconducting 
is  given  by 

H  =  am  +  ECN2  -  JS{S  +  1)  +  ggbSB.  (7) 

i 

Note  that  the  only  random  component  of  the  problem 
is  the  single-particle  spectrum  £j,  while  the  exchange  J 
and  charging  energy  Ec  do  not  fluctuate.  It  should  be 
realized,  though,  that  Eqs.  (H|7j)  become  exact  only  in 
the  limit  gr  — >  oo.  The  corrections  to  Eqs.  ([|ff|),  which 
appear  at  finite  gr  are  sometimes  important.  However, 
if  <7t  1,  these  corrections  do  not  bring  essentially  new 

physics  to  the  problem  of  spin  magnetization. 

The  conductance  is  calculated  using  the  many-particle 
energies  Ejy  and  wave  functions  obtained  numerically  for 
GOE  and  GUE  random  matrix  realizations.  An  example 
of  the  peak  positions  and  peak  spacing  evolution  as  func¬ 
tions  of  the  magnetic  field  for  a  particular  GOE  realiza¬ 
tion  is  presented  in  Fig.  ^]for  (a)  J  =  0  and  (b)  J  =  O.lJi. 
For  both  cases  we  present  g'N  =  —  Ec.  The  noninter¬ 

acting  case,  Fig.  [tja,  shows  all  the  previously  described 
features  of  the  Pauli  behavior.  Once  a  weak  exchange  in¬ 
teraction  is  included  the  behavior  changes  qualitatively, 
not  just  at  high  field,  but  also  in  the  vicinity  of  B  =  0.  In 
particular,  peak  positions  are  not  always  paired,  result¬ 
ing  in  occasional  field-independent  peak  spacings.  This 
occurs  when  two  consecutive  orbitals  are  first  filled  with 
down  spin  electrons  and  only  later  they  acquire  up  elec¬ 
trons.  Generally  the  enhancement  of  the  spin  of  the  dot 
by  S  will  be  accompanied  by  a  2 S  bunch  peaks  moving 
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with  the  same  slope.  If  the  peak  spacing  is  plotted, 
two  sets  of  25  —  1  flat  curves  sandwiching  a  sloped  one 
will  appear.  Changes  to  the  crossing  also  attributed  to 
a  sudden  change  of  the  GS  spin  which  is  not  associated 
with  a  single  electron  orbital  crossing. 

As  illustrated  above,  the  appearance  of  higher  spin 
clearly  manifests  itself  in  the  peak  position  and  peak 
spacing  trajectories.  Using  Eq.  (f?)  one  can  predict  the  fre¬ 
quency  of  spontaneous  magnetization  appearances.  For 
weak  exchange,  J  -C  Si ,  the  probability  that  the  dot  GS 
has  a  spin  S  is  determined  by  the  probability  of  finding 
2 S  orbitals  so  close  to  each  other  that  the  gain  in  ex¬ 
change  energy  (due  to  the  polarization)  overwhelms  the 
loss  in  the  kinetic  energy  (due  to  single  rather  than  dou¬ 
ble  occupations  of  the  orbitals). 

Thus  the  probability  for  a  certain  value  of  ground  state 
magnetization  S  boils  down  to  the  question  of  the  prob¬ 
ability  of  finding  2 S  close  single  electron  levels.  That 
probability  may  be  estimated  using  random  matrix  the¬ 
ory,  as  we  now  present. 

The  probability  of  finding  a  set  of  single  electron  or¬ 
bital  at  energies  is  given  by  Eq.  (||).  Thus  as  long 
as 


£i+2S  —  £i  <  JS(S  +  1)  +  gfiBBS,  (8) 


a  ground  state  spin  of  at  least  S  will  appear.  Taking  into 
account  the  probability  of  Eq.  (|j),  it  is  possible  to  write 
the  probability  of  obtaining  a  ground  state  spin  S'  as  a 
function  of  the  linear  combination: 


y  _  J  .  9VbB 

6 i  |S+  iJV 


(9) 


resulting  in  the  following  probability 


One  may  wonder  how  well  does  the  above  picture  hold 
for  systems  for  which  real  correlation  between  the  elec¬ 
tron  exist,  and  the  dimensionless  conductance  is  not  too 
large.  A  canonical  example  for  such  systems  is  the  Hub¬ 
bard  model: 

H  =  eiai,<Tai,cr  —  t  yy  (ai,crai,o-  +  h.c) 
i,a  <i,j>,a 

+U  ^2  aI,tabfai,4,ab!  +  Q^bSB  (11) 


where  <  i,j  >  denotes  nearest  neighbor  lattice  site,  a\  a 
is  an  creation  operator  of  an  electron  at  site  i  with  spin  er, 
ei  is  the  site  energy,  chosen  randomly  between  —W/2  and 
W/2  with  uniform  probability,  and  U  is  the  interaction 
constant. 

The  Hilbert  space  is  large  even  for  relatively  small  sys¬ 
tems.  For  example  for  a  4  x  4  lattice  with  6  electrons 
for  the  S  =  0  sector  has  the  size  of  313600.  Using  the 
Lanczos  method  we  obtain  the  many-particle  eigenval¬ 
ues  En(S)  as  function  of  the  ground  state  in  each  spin 
sector  S  for  different  values  of  interaction  U.  The  dis¬ 
order  was  chosen  as  W  =  8 1,  which  corresponds  to  the 
metallic  regime,  although  the  value  of  the  dimensionless 
conductance  is  quite  low  gx  ~  0(1). 

We  calculated  the  probability  for  the  appearance  of  a 
specific  value  of  GS  spin  S  for  different  values  of  U  and 
magnetic  field  B.  In  order  to  use  the  scaling  parameter 
X  we  need  to  deduce  the  appropriate  value  of  J  for  each 
value  of  U.  This  has  been  done  by  fitting  the  dependence 
of  the  average  lowest  energy  in  a  given  spin  sector  {E{S)) 
to 


{E(S)}  =  Sis2  -  JS(S  +  1).  (12) 


Pj,b(S)  =  Cfx(/3S+1)(2S'-1)(l  -  K%X2),  (10) 

The  coefficients  Cg  and  Kg  depend  on  both  (3  and  S. 
Their  values  for  S  =  1,3/2  are  presented  in  table  I. 


(3  =  1 

(3  =  2  I 

S=  1 

S  =  3/2 

5  =  1 

5  =  3/2 

c 

ti-73 

9tt4/50 

87/9 

817/400 

K 

7T2/5 

187/49 

87/25 

7927/1225 

TABLE  I.  The  factors  Cg  and  Kg  appearing  in  Eq.  ([loj). 


Thus  the  RMT  model  presents  a  striking  conclusion. 
The  influence  of  exchange  interactions  and  externally  ap¬ 
plied  magnetic  field  on  the  ground  state  magnetization  of 
a  quantum  dot  may  be  summed  up  in  a  single  parameter 
scaling  function,  where  the  scaling  parameter  is  a  simple 
linear  combination  of  J  and  B.  As  has  been  shown  from 
numerical  simulations  |23),  the  scaling  holds  for  larger 
values  of  X  than  expected  from  the  first  and  second  term 
perturbative  analysis  presented  in  Eq.  ©• 


The  distribution  Pj(u),b(S)  based  on  numerical  results 
for  1000  different  realizations  of  the  Hubbard  model,  for 
different  values  of  U  and  B  for  (3  =  1  and  (3  =  2  is  pre¬ 
sented  in  Fig.  |}  Although  some  deviations  do  appear, 
especially  for  the  higher  values  of  U ,  one  can  see  that 
the  overall  form  of  the  scaling  function  holds  remarkably 
well.  Thus,  though  the  Hamiltonian  in  Eq.  (iT)  is  obtained 
under  the  conditions  of  large  gx  and  no  correlations  in 
the  system,  it  nonetheless  appears  to  capture  much  of 
the  physics  even  for  a  moderate  gx  and  in  the  presence 
of  electron-electron  correlations.  The  reason  for  this  is 
that  although  neither  J  nor  Ec  is  constant  once  gx  ~  1, 
their  fluctuations  remain  small  relative  to  their  average, 
as  is  demonstrated  in  Fig.  |j. 


III.  EXPERIMENT 

In  this  section  we  describe  recent  experimental  mea¬ 
surements  of  GS  spin  for  a  gate-defined  GaAs  quantum 
dot  containing  roughly  400  electrons.  As  discussed  in  the 
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X 

FIG.  2.  Probabilities  of  S  =  1  for  J  =  O.ldi  as  functions  of 
the  magnetic  field  B  for  (a)  the  GOE  and  (b)  GUE  cases.  The 
small  symbols  represent  numerical  results  for  the  Hubbard 
model,  the  large  filled  symbols  represent  numerical  results  for 
the  RMT  model,  while  the  curves  (solid  -  first  order,  dashed 
-  second  order)  represent  Eq.  (|Ioj). 


Introduction,  the  dot  is  coupled  to  electron  reservoirs  via 
tunnelling  leads  (i.e. ,  g  <  2 e2//i  for  both  leads)  so  trans¬ 
port  is  dominated  by  CB  effects.  Measurements  were 
carried  out  at  sufficiently  low  temperature  and  bias  that 
the  differences  in  GS  energies  between  dots  with  N  +  1 
and  N  electrons  should  be  extractable  from  CB  peak  po¬ 
sitions,  VgN\ 

The  dot  is  formed  at  the  interface  of  a  GaAs/AlGaAs 
heterostructure  (90  nm  below  the  wafer  surface)  by  elec¬ 
trostatic  depletion  using  surface  gates.  A  Si  delta  dop¬ 
ing  region  is  located  40  nm  above  the  heterointerface. 
The  two-dimensional  electron  gas  (2DEG)  has  density 
~  2.0  x  1011  cm-2  and  bulk  mobility  ~  1.4  x  105  cm2 /Vs, 
yielding  a  transport  mean  free  path  ~  1.5  /im.  The 
small  dot  area,  A  ~  0.25  /im2,  makes  transport  predom¬ 
inantly  ballistic  within  the  device.  Characteristic  en¬ 
ergy  scales  for  the  measured  device  include  the  mean 
level  spacing,  A  =  2nh2/m*A  ~  30  g.eV,  the  charg¬ 
ing  energy,  Ec  ~  400/reV,  and  the  Thouless  energy 
Eth  =  TivfA~ 1/2  ~  340/reV.  Measurements  were  carried 
out  in  a  dilution  refrigerator  with  a  mixing  chamber  tem¬ 
perature  of  25  mK  using  standard  ac  lock-in  techniques 
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FIG.  3.  The  exchange  parameter  J  as  functions  U  for  a 
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with  a  source-drain  bias  voltage  of  2  fieV .  A  base  elec¬ 
tron  temperature  of  Te  ~  50  mK  was  determined  from 
CB  peak  widths. 

To  allow  the  magnetic  field  to  couple  predominantly  to 
spin,  the  sample  was  oriented  with  the  plane  of  the  elec¬ 
tron  gas  along  the  axis  of  the  primary  solenoid,  aligned 
manually  to  within  0.5  degrees.  In  addition,  a  small  pair 
of  coils  attached  to  the  vacuum  can  of  the  fridge,  ori¬ 
ented  perpendicular  to  the  plane  of  the  sample  (see  Fig. 
4),  was  used  to  null  out  any  perpendicular  held  from 
misalignment  as  well  as  to  explicitly  break  time-reversal 
symmetry.  Both  primary  and  trimming  coils  were  un¬ 
der  computer  control,  allowing  sweeps  of  strictly  paral¬ 
lel  held  by  trimming  out  any  perpendicular  component 
for  any  applied  parallel  held.  (We  estimate  the  uncer¬ 
tainty  in  B±  to  be  less  than  </>0/4  through  the  dot  at 
B\\  =  5  T.)  Despite  this  precise  held  trimming  capabil¬ 
ity,  similar  measurements  j(|  in  larger  dots,  fabricated 
on  the  same  wafer,  indicate  orbital  coupling  due  to  the 
strictly  parallel  held  for  B\\  >  OAT  (evident,  for  instance, 
from  the  disappearance  of  the  weak  localization  feature 
at  B±  =  0  at  higher  B\\).  The  origin  of  this  surprising 
coupling — which  makes  interpretation  of  our  high-held 
data  difficult — is  still  under  investigation. 

Conductance  measurements  across  ten  consecutive 
Coulomb  blockade  peaks,  measured  as  a  function  of  Vg 
and  B\\  (i.e.,  strictly  By,  properly  trimmed),  are  shown  in 
Fig.  |a.  More  positive  gate  voltage  corresponds  to  higher 
energy,  and  can  be  calibrated  from  the  CB  “diamonds” 
using  high  source-drain  bias  measurements  [|l).  All  data 
are  taken  with  B±  =  20 mT  in  order  to  ensure  that  time- 
reversal  symmetry  is  broken,  which  changes  the  statistics 
of  dot  wave  functions. 

In  addition  to  the  individual  motions  of  each  peak, 
there  is  a  diamagnetic  shift  of  all  peaks  (see  Fig.  &>). 
presumably  due  to  the  effect  of  the  parallel  held  on  the 
effective  well  conhnement  potential.  The  slight  param¬ 
agnetic  shift  visible  at  low  held,  (J3||  <  OAT),  is  not 
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(a)  (b)  _^-4 


FIG.  4.  a)  Schematic  diagram  of  a  gate-defined  quantum 
dot,  indicating  the  orientations  of  B y  and  B±  relative  to  the 
planar  quantum  dot.  (Orientation  of  B y  within  the  plane  is 
not  accurately  depicted  and  has  not  been  investigated.)  b) 
Diagram  showing  placement  of  superconducting  coils  used  to 
generate  B±  on  vacuum  can,  inside  primary  solenoid  used  to 
generate  Bn. 

understood  at  present.  In  the  following  analysis  of  data, 
the  common  curve  (in  Fig.  |b)  has  been  subtracted  from 
each  peak  position. 

Peak  positions  and  peak  spacings  extracted  from  the 
peaks  in  Fig.  5  are  shown  in  Fig.  6.  The  slope  of  peak 
positions  as  a  function  of  B ||  is  consistent  with  a  Zeeman 
energy  term  Es  =  ±^ggsB,  using  the  ^-factor  for  bulk 
GaAs,  |g|  =  0.44.  As  discussed  in  the  Introduction,  alter¬ 
nating  slopes  for  consecutive  peaks  would  indicate  an  al¬ 
ternating  0,  0,  5...  GS  spin  structure.  Our  data,  on  the 

other  hand,  shows  three  consecutive  pairs  of  peaks  mov¬ 
ing  with  the  same  slope,  suggesting  the  presence  of  higher 
spin  states.  Proposed  values  for  the  eight  consecutive  GS 
spin  states  shown  here  are  included  in  Fig.  6a.  We  em¬ 
phasize,  however,  that  these  are  only  plausible  values  for 
the  spin;  it  is  not  possible  to  determine  unambiguously 
the  absolute  magnitude  of  GS  spin  from  measurements 
of  peak  position,  which  reflect  changes  in  spin  from  the 
N  to  N+l  ground  states.  The  values  shown  in  Fig.  6a 
minimize  the  ground  state  spin  for  the  system. 

In  the  proposed  spin  labelling  scheme,  three  out  of  the 
five  even-iV  states  have  5=1,  i.e. ,  P{S  =  1)  ~  0.6. 
This  fraction  of  5  =  1  relative  to  the  number  of  5  =  0 
states  is  well  beyond  the  expected  value  given  reasonable 
estimates  of  J  for  GaAs  dots.  Of  course,  with  only  five 
spin  states  considered,  statistics  are  quite  poor.  Further 
experiments  are  needed  to  see  if  this  discrepancy  is  sig¬ 
nificant. 

Figure  7  shows  that  peak  spacings  clearly  separate  into 
three  branches,  a  top  branch  with  slope  roughly  gg  (cor¬ 
responding  to  a  GS  spin  decrement  followed  by  an  in¬ 
crement)  a  bottom  branch  with  slope  roughly  —gg  (cor¬ 
responding  to  a  GS  spin  increment  followed  by  a  decre¬ 
ment)  and  a  middle  branch  with  slope  near  zero  (corre¬ 
sponding  to  two  consecutive  increments  or  decrements). 
The  existence  of  the  middle  branch  is  the  signature  of 
higher  GS  spins.  The  good  agreement  between  the  slopes 
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FIG.  5.  a)  Nine  consecutive  Coulomb  blockade  peaks  mea¬ 
sured  as  a  function  of  gate  voltage  Vg  and  parallel  field  B. 
Conductance  is  shown  in  grayscale,  with  black  indicating  high 
conductance  and  white  low  conductance.  The  field  shown  on 
the  bottom  axis  is  the  field  strictly  parallel  to  the  heterointer¬ 
face;  perpendicular  field  is  held  constant  at  20 mT  as  discussed 
in  the  text.  Note  that  conductance  in  the  valleys  is  low  but 
not  zero,  due  to  strong  tunnelling  in  the  leads,  b)  Average 
peak  position  of  the  peaks  in  a),  showing  a  diamagnetic  shift 
to  higher  energies  common  to  all  peaks.  This  average  is  sub¬ 
tracted  from  all  peak  positions  before  further  analysis. 

of  the  upper  and  lower  branches  and  the  expected  slopes 
of  gg,  as  well  as  the  absence  of  a  range  of  intermedi¬ 
ate  slopes  suggest  that  the  peak  spacing  reflects  spin 
rather  than  orbital  coupling.  At  higher  fields,  the  di¬ 
rections  of  peak  motion  change,  often  abruptly  and  from 
one  straight  segment  to  another,  as  seen  in  Fig.  8.  This 
behavior  is  qualitatively  similar  to  the  numerical  data  in 
Fig.  lb.  The  rounding  of  straight  segments  where  the 
slope  changes  presumably  results  from  spin-orbit  inter¬ 
action  which  mixes  spins,  and  may  provide  a  direct  mea¬ 
sure  of  spin-orbit  interactions  in  dots.  It  is  interesting  to 
note  that  the  peaks  analyzed  in  Figs.  6  through  8  were 
measured  in  a  regime  of  high  tunneling  conduction  in  the 
leads.  This  can  be  seen  by  noting  the  grayscale  of  Fig.  5. 
When  the  dot  is  more  pinched  off  from  the  reservoirs,  so 
that  the  CB  peaks  have  a  height  of  0.1  e2 /h  or  less,  peak 
motion  is  more  difficult  to  interpret,  and  does  not  seem 
to  follow  the  clear  patterns  illustrated,  for  instance,  in 
Fig.  8.  As  seen  in  Fig.  9,  peak  spacings  for  four  consecu¬ 
tive  weak-tunneling  peaks  at  do  not  move  with  slopes  of 
gg  at  low  B\\ .  Though  the  data  are  noisier  in  this  con¬ 
figuration,  it  is  clear  that  the  slopes  are  in  all  cases  less 
than  gg.  We  do  not  have  an  explanation  for  the  different 
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FIG.  6.  a)  Peak  positions  (with  average  motion  subtracted) 
as  a  function  of  parallel  field  for  the  nine  peaks  shown  in  Fig. 
5a.  Straight  lines  indicate  expected  peak  motion  for  spin-i 
transitions  using  the  g  factor  for  bulk  GaAs,  |p|  =  0.44.  Pos¬ 
itive  slopes  indicate  transitions  to  lower  spin  states;  negative 
slopes  indicate  transitions  to  higher  spin  states.  Numbers 
indicate  a  possible  ground-state  spin  structure  for  ten  consec¬ 
utive  ground  states.  Peak  position  data  are  offset  for  clarity, 
b)  Peak  spacings  for  the  position  data  shown  in  a),  offset  for 
clarity.  Solid  lines  indicate  expected  motion  of  peak  spacing 
for  spin--);  transitions,  using  |g|  =  0.44 

behavior  depending  on  if  the  dot  is  more  closed  off  or 
nearly  open. 

IV.  DISCUSSION 

In  conclusion,  we  have  shown  that  a  relatively  weak  ex¬ 
change  interaction  qualitatively  explains  the  deviations 
from  the  Pauli  picture  seen  in  recent  experiments  ||[f],[| . 
We  calculate  the  probability  that  different  values  of  spin 
appear  as  a  function  of  exchange  interaction  and  mag¬ 
netic  field  coupling  to  the  spin.  The  strength  of  the 
exchange  interaction,  J,  and  the  effective  g  factor  are 
the  only  adjustable  parameters  that  determine  the  prob¬ 
abilities  for  a  GS  of  the  dot  to  have  any  given  spin  at 
any  given  magnetic  field  for  both  GOE  and  GUE  cases. 
In  particular  we  predict  that  these  probabilities  will  fol¬ 
low  a  one-parameter  scaling  law  over  a  wide  range  of 
magnetic  fields  and  exchange  interaction  strengths.  We 


FIG.  7.  Peak  spacings  from  Fig.  6b,  offset  to  align  spacings 
at  B ||  =  0.  Slopes  cluster  in  three  branches,  interpreted  as  a 
change  from  decreasing  to  increasing  spin  transitions  (positive 
slope),  two  consecutive  spin  transitions  in  the  same  direction 
(zero  slope),  and  a  change  from  increasing  to  decreasing  spin 
transitions  (negative  slope).  Solid  lines  indicate  expected  mo¬ 
tion  of  peak  spacing  for  spin-i  transitions,  using  |gr|  =  0.44, 
with  no  adjustable  parameters. 

also  report  preliminary  experiments  to  investigate  the  ex¬ 
change  effects  in  the  ground  state  spin  of  quantum  dots. 
The  signature  of  interaction  effects  is  the  appearance  of 
higher-spin  ground  states,  which  show  up  in  the  experi¬ 
mental  data  as  peak  spacing  traces  that  have  zero  slope 
as  a  function  of  parallel  magnetic  field.  Such  features  are 
indeed  observed,  as  seen  for  instance  in  Fig.  8.  Further 
experiments  are  needed  to  obtain  sufficient  statistics  to 
quantitatively  test  the  predictions  of  theory. 
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